ASTRONOMY
AND
ASTROPHYSICS
February 1, 2008
A&A manuscript no.
(will be inserted by hand later)
Your thesaurus codes are:
02.12.3 ;11.09.1 Mkn 110; 11.19.1
Strong optical line variability in Mkn 110
a\
0\
K. Bischoff and W. Kollatschny
Universitats-Sternwarte Gottingen, Geismarlandstrafie 11, D-37083 Gottingen, Germany
Received date; accepted date
>
(N
o
CO
O
as
a\
Oh'
I 1
O 1
X
Abstract. We present results of a long-term variability
campaign on the Seyfert 1 galaxy Markarian 110. Mkn 110
is a narrow-line Seyfert 1 object hosted in a morphologi-
cal peculiar galaxy. We monitored the optical continuum
and the line intensities as well as their profiles over a
time interval of nearly ten years. The continuum and the
Balmer lines varied by a factor of 2 to 5 within two years.
The HeIIA4686 line showed exceptional intensity varia-
tions of a factor of eight. We detected an additional inde-
pendent very broad-line region in high intensity stages of
the Balmer and Hell lines. The CCF analysis of the Hell
line indicates that this very broad-line region originates at
a distance of 9 light days only from the central ionizing
source.
Key words: galaxies: individual:Mkn 110 - galaxies: Sey-
fert - lines: profiles
1. Introduction
Markarian 110 is a nearby (z=0.0355) Seyfert 1 galaxy
with highly irregular morphology. The apparent magni-
tude of the total system is my=15.4 mag (Weedman 1973)
corresponding to My=-20.4 mag (Hq—75 km s _1 Mpc -1 ).
A foreground star is projected on the host galaxy at a dis-
tance of 6 arcsec to the nucleus in north-east direction (see
Fig. 1). Therefore, it has been supposed in some early
papers that Mkn 110 might be a double nucleus galaxy
(Petrosian et al. 1978).
On the other hand the peculiar morphology of Mkn 110
is an indication for a recent interaction and/or merging
event in this galaxy (Hutchings & Craven 1988). One can
clearly recognize a tidal arm to the west with a projected
length of 50 arcsec (corresponding to 35 kpc) and further
Send offprint requests to: W. Kollatschny
* Visiting Astronomer, German-Spanish Astronomical Cen-
tre, Calar Alto, operated by the Max-Planck-Institute for As-
tronomy, Heidelberg, jointly with the Spanish National Com-
mission for Astronomy.
Correspondence to: kbischo@gwdg.de
o ©
■
«
j i I i i i I i i i I i i i I i i
-SO 30 40
RA [arcsec]
Fig. 1. Optical R-band image of Mkn 110. North is up,
east is to the left.
signs of asymmetry in the disturbed host galaxy on the
R-band image (Fig. 1).
Ten years ago we started a long-term variability cam-
paign to study the continuum and emission line intensity
variations in selected AGNs. Besides our principal inter-
est in the long-term variability behaviour of these galaxies
on its own we want to compare the individual variability
properties with those of other galaxies from the interna-
tional AGN watch campaign (Peterson et al. 1991) (e.g.
NGC 5548, Kollatschny & Dietrich 1996) and LAG cam-
paign (Robinson 1994) (e.g. NGC 4593, Kollatschny &
Dietrich 1997). Further results on continuum and H/3 vari-
ations in Mkn 110 have been published in a recent paper
on variability of Seyfert 1 galaxies (Peterson et al. 1998a).
2. Observations and data reduction
We took optical spectra of Mkn 110 at 24 epochs from
February 1987 until June 1995. The sampling of the ob-
2
K. Bischoff & W. Kollatschny: Strong optical line variability in Mkn 110
servations extends from days to years. In Table 1 we list
our observing dates and the corresponding Julian Dates.
The spectra were obtained at Calar Alto Observatory in
Spain with the 2.2 m and 3.5 m telescopes as well as at
McDonald Observatory in Texas with the 2.1m and 2.7 m
telescopes. The individual exposure times range from 10
minutes to 1 hour (see Table 1). We used spectrograph
slits with projected widths of 2 to 2.5 arcsec and 2 arcmin
length under typical seeing conditions of 1 to 2 arcsec.
We extracted spectra of the central 5 arcsec. The slit was
oriented in north-south direction in most cases.
To investigate the spatial extension of the narrow line
region we took spectra at different position angles: 0°, 45°,
90°, and 135°. A possible extended [OIII]A5007 emission
line flux was always less than 3% of the nuclear point-like
emission. Mkn 110 has been inspected by Nelson et al.
(1996) with the Hubble Space Telescope WFPC-1 in the
near infrared spectral range. They detected a dominant
unresolved nucleus in this galaxy.
Our optical spectra typically cover a wavelength range
from 4000 A to 7200 A with a spectral resolution of 3 to
7A per pixel. We used different CCD detectors in the
course of this monitoring program: until 1989 a RCA-
chip (1024x640), in January and July 1992 a GEC-
chip (1155x768), and in August 1992 a Tektronix-chip
(1024x1024).
The reduction of the spectra (flat fielding, wavelength
calibration, night sky subtraction, flux calibration, etc.)
was done in a homogeneous way using the ESO MIDAS
package.
The absolute calibration of our spectra was achieved by
scaling the [OIII]A5007 line of all spectra to those obtained
under photometric conditions. Our absolute [OIII]A5007
flux corresponds within 5% to that obtained by Peter-
son et al. (1998a). For a better comparison of these two
data sets we will use exactly the same [OIIIJA5007 flux
of 2.26 10~ 13 ergs -1 cm~ 2 . Furthermore, we corrected all
our data for small spectral shifts and resolution differ-
ences with respect to a mean reference spectrum using an
automatic scaling program of van Groningen & Wanders
(1992).
The R-band image of Mkn 110 (Fig. 1) was taken with
the 2.2m telescope at Calar Alto observatory on Septem-
ber 20, 1993 with an exposure time of 6 minutes. Again, we
reduced this CCD image with the ESO MIDAS package.
In the course of our discussion we additionally will
make use of archival IUE spectra of Mkn 110 taken on
February 28 and 29, 1988.
3. Results
Some typical spectra are plotted in Fig. 2 showing the
range of intensity variations. Immediately one can recog-
nize the strong variations in the continuum, in the Balmer
lines and especially in the HeIIA4686 line. The continuum
variations are most pronounced in the blue section. The
Table 1. Log of observations
Julian Date
U l uate
Telescope
Exp. time
Z 44U UUU+
[sec]
( 2 )
(A\
( 4 )
6828
1987-02-01
CA 3.5
1500
7229
1988-03-08
CA 3.5
4800
7438
1988-10-03
CA 3.5
1200
7574
1989-02-16
CA 2.2
3000
7663
1989-05-16
CA 3.5
1800
7828
1989-10-28
CA 2.2
1200
8632
1992-01-09
CA 3.5
600
8812
1992-07-08
CA 3.5
1800
8860
1992-08-25
CA 2.2
900
8862
1992-08-27
CA 2.2
900
8864
1992-08-25
CA 2.2
900
9078
1993-03-31
CA 2.2
3600
9080
1993-04-02
CA 2.2
3600
9083
1993-04-05
CA 2.2
2400
9123
1993-05-15
MDO 2.1
1800
9237
1993-09-06
CA 3.5
1200
9419
1994-03-07
MDO 2.1
2100
9595
1994-08-30
CA 2.2
1200
9776
1995-02-27
MDO 2.1
1800
9785
1995-03-08
MDO 2.1
1200
9786
1995-03-09
MDO 2.1
1800
9787
1995-03-10
MDO 2.1
1200
9811
1995-04-03
MDO 2.1
1200
9870
1995-06-01
MDO 2.7
900
CA = Calar Alto Observatory
MDO = McDonald Observatory
4000 5000
a [A]
Fig. 2. Normalized spectra of Mkn 110 taken at different
epochs in Oct. 1988, Oct. 1989, March 1988, Feb. 1989,
May 1989, Jan. 1992 (from bottom to top).
continuum gradient changes as a function of intensity. The
emission line profiles of the Balmer lines in Mkn 110 are
quite narrow (FWHM(H/3)=1800 km s" 1 ) similar to those
of the so called narrow-line Seyfert 1 galaxies. Week Fell
emission is present in the spectra blending the red line
K. Bischoff & W. Kollatschny: Strong optical line variability in Mkn 110
3
4000
5000
^ [A]
6000
7000
Fig. 3. Difference spectra with respect to our minimum
stage in Oct. 1988; otherwise same epochs as in Fig. 2.
wings of [OIII]A5007 and H/3. We measured the integrated
intensity of Fell line blends between 5134 A and 5215 A.
The main Fell components in this region are the 5169 A
and 5198 A lines belonging to the multiplets 42 and 49.
The Fell line flux (2.0 10" 14 ergs
) remained con-
stant during our variability campaign within the error of
10%.
Difference spectra with respect to our minimum stage
in October 1988 are plotted in Fig. 3. All narrow line
components cancel out. The Fell lines disappear in the
difference spectra as well. In the Balmer profiles (e.g. H/3)
very broad, slightly rcdshifted components stand out in
the high intensity stages.
3.1. Line and continuum variations
The results of our continuum intensity measurements at
3750 A, 4265 A, and 5100 A as well as the integrated line
intensities of Ha, H/3, HeIIA4686, HeIA5876, and HcIA4471
are given in Table 2. The individual light curves are plot-
ted in Fig. 4.
The continuum intensities are mean values of the wave-
length ranges given in Table 3, column (2). Line intensities
were integrated in the listed limits after subtraction of a
linear pseudo-continuum defined by the boundaries given
in column (3). All wavelengths arc given in the rest frame.
We started our monitoring program in 1987. There-
fore, our 5100 A continuum light curve covers up a larger
time interval than the light curve of Peterson et al. (1998a)
beginning in 1992. The observing epochs are partly com-
plementary in the common monitoring interval. But, spec-
tra taken nearly simultaneously from both groups (within
one week) correspond with each other to better than 5%
in the continuum fluxes. In Fig. 5 we compare our con-
tinuum light curve with that of Peterson et al. (1988a) in
the common interval of observations. Both light curves are
Table 3. Boundaries of mean continuum values and line
integration limits
Cont./Line Wavelength range Pseudo-continuum
(1) (2) (3)
Cont. 3750
3745 A
- 3755 A
Cont. 4265
4260 A
- 4270 A
Cont. 5100
5095 A
- 5105 A
HeIIA4686
4600 A
- 4790 A
4600 A
- 4790 A
HeIA4471
4430 A
- 4530 A
4265 A
- 4600 A
HeIA5876
5800 A
- 5960 A
4600 A
- 5100 A
H/3
4790 A
- 4935 A
4600 A
- 5100 A
Ha
6420 A
- 6770 A
6420 A
- 6770 A
48750 49000 49250 49500 49750
JD (+2400000)
Fig. 5. Comparison of our continuum light curve at 5100 A
(filled circles) with that of Peterson et al. (1988a) between
JD 2448600 and JD 2500000.
in very good accordance regarding to the intensity varia-
tions. Fig. 5 shows that the data sets are not significantly
undersampled. However, our H/3 intensities are system-
atically higher than those of Peterson et al. (1998a) as
we integrated over a larger wavelength range and we car-
ried out a slightly different continuum subtraction. This
method led to a lower pseudo-continuum flux at 4790 A.
The H/3 fluxes are in perfect agreement if we multiply the
values given by Peterson et al. (1998a) by a factor of 1.15.
The pattern of the continuum light curves at 5100 A
and 4265 A (Fig. 4) is identical apart from their differ-
ent amplitudes. The HeIIA4686 light curve follows closely
these continuum light curves. The light curves of the
Balmer lines Ha and H/3 are similar among themselves
and the light curves of both Hcl lines as well.
Statistics of our measured continuum and emission line
variations are presented in Table 4. We list minimum and
maximum fluxes F min and F max , peak-to-peak amplitudes
F^max = F max /F min , the mean flux over the entire period
4
K. Bischoff & W. Kollatschny: Strong optical line variability in Mkn 110
i
i
48000 49000
JD (+2400000)
47000 48000 49000
JD (+2400000)
250
.1 :
200
150
I
J :
if i :
100
Y A , /> :
50
'i :
48000
JD (+2400000)
! s
48000
JD (+2400000)
'i f
48000 49000
JD (+2400000)
G
O
- 500
i i i i
i
; I.
i
i
i,-
f
''' :
''i
¥
i i i i
i
48000
JD (+2400000)
= 10 -
48000 49000
JD (+2400000)
Fig. 4. Light curves of continuum flux at 5100 A and 4265 A (in units of 10 15 erg cm 2 s 1 A x ) and integrated
emission line flux of Ha, H/3, HeIIA4686, HeIA5876 and HeIA4471 (in units of 10" 15 erg cm" 2 s" 1 ). The points arc
connected by a dotted line to aid the eye.
K. Bischoff & W. Kollatschny: Strong optical line variability in Mkn 110
Table 4. Variability statistics
Cont./Line
(1)
Fmm
(2)
(3)
(4)
<F>
(5)
o F
(7)
Cont. 3750
2.94
15.10
5.14
7.51
4.322
0.573
Cont. 4265
1.96
7.56
3.86
3.91
1.533
0.391
Cont. 5100
1.53
5.93
3.88
3.26
1.044
0.318
HeIIA4686
30.6
237.6
7.76
102.6
48.95
0.466
HeIA4471
12.0
41.7
3.47
28.56
7.19
0.180
HeIA5876
33.1
103.2
3.12
69.44
17.31
0.236
H/J
266
624.5
2.34
419.2
79.94
0.189
Ha
876
2091
2.39
1572
265.3
0.168
Continuum flux in units of 10
-15 —1 -
erg sec cm
2 A- 1 .
Line flux in units of 10 15 erg sec 1 cm
of observations <F>, the standard deviation a F , and the
fractional variation
< F >
as defined by Rodriguez-Pascual et al. (1997). The ex-
treme variability amplitudes of Mkn 110 attract attention
compared to other galaxies, e.g. NGC 4593 (Dietrich, Kol-
latschny et al. 1994) and the Seyfert 1 galaxies from the
sample of Peterson et al. (1998a).
The variability amplitudes of the continuum increase
towards the short wavelength region. These amplitudes as
well as those of the emission line intensities are exception-
ally high. The variability amplitude of the HellA4686 line
is unique compared to the other emission lines in Mkn 110
and compared to optical lines in other Seyfert galaxies.
In Fig. 6 we plot the line intensity ratios HeIIA4686/H/3
and HeIA5876/H/3 as a function of continuum intensity at
5100 A. These line intensity ratios increase slightly for the
Hcl line but strongly for the highly ionized Hell line.
3.2. Balmer decrement
We calculated Balmer decrement Ha/H/3 values in the
range from 3.2 to 4.3. Simple photoionization calculations
(Case B) result in a value of 2.8 for this line ratio (Oster-
brock 1989). Deviations of the observed Balmer decrement
from the theoretical value are often explained by wave-
length dependent dust absorption and/or by collisional
excitation effects. We will show later on that the observed
difference can not be explained by dust absorption alone
in the broad- line region clouds of Mkn 110. There is a clear
anti-correlation of the Balmer decrement with the contin-
uum flux (Fig. 7). One has to keep in mind that the indi-
vidual Balmer lines and the continuum of each spectrum
originate in distinct regions of the BLR at different times.
This will be confirmed by the cross-correlation analysis
later on. A very tight correlation can not be expected be-
cause of the short-term variations in this galaxy. The solid
line in Fig. 7 is a linear fit to all our data points except
0.4
0.3
x
ac
u. 0.2
o.i
+
+
r
i
2 3 4 5 6
F 5100 [10" 15 erg s" 1 cm" 2 A" 1 ]
Fig. 6. Line intensity ratios HeIIA4686/H/3 (filled circles)
and HeIA5876/H/3 (open circle) as a function of continuum
flux at 5100 A.
4.4
3G
fc. 3.8
a
3G
b. 3.6
2 3 4 5 6
F 5100 [10" 15 erg s" 1 cm" 2 A" 1 ]
Fig. 7. The Balmer decrement Ha/H/3 as a function of
continuum flux at 5100 A. The minimum value of F 5100
was omitted for the linear fit.
for the lowest continuum intensity point. At this epoch
(JD +5959) the Ha intensity was extreme low (Fig. 4) in
contrast to H/3 and the continuum. An anti-correlation of
the Balmer decrement with the continuum flux has been
first noted in NGC 4151 by Antonucci & Cohen (1983).
3.3. UV spectra
Two UV spectra haven been taken with the IUE satellite
with a time interval of one day only. These spectra have
been taken nearly simultaneously (within 8 days) to our
optical observations in March 1988. Therefore, these two
spectra are suitable for a determination of optical/UV line
400
300
to 200
u
<u
S 100
p2=
1200
K. Bischoff & W. Kollatschny: Strong optical line variability in Mkn 110
l
x [A]
1600
1800
0.8
0.6
h 0.4
0.2
-0.2
Cont5100 x
-200
-100 100
t [days]
200
Fig. 8. Short
taken on Feb.
wavelength IUE UV spectra of Mkn 110 Fig. 9. Cross-correlation functions CCF(r) of HeIIA4686,
28 and Feb. 29, 1988
HeIA5876, H/3 and Ha light curves with the 5100 A con-
tinuum light curve.
intensity ratios. Fig. 8 shows an overplot of both short
wavelength UV spectra taken with an time interval of 1
day. The spectra are identical in the continuum and in
the emission lines within the error limits. The different
flux values in the center of the CIVA1550 line are due to
saturation effects in one of the spectra.
We determined an integrated Lya flux of (4.2 ±
0.3) 10~ 12 ergs -1 cm~ 2 . Comparison with the optical spec-
tra results in a Lya/H/3 ratio of 11.0 at the observing
epoch March, 1988.
The HcIIA1640 flux amounts to (1.4 ± 0.2)10" 13
ergs" 1 cm" 2 . The HeIIA1640/HeIIA4686 ratio of 2.45 is
about a factor of two lower than that of typical photoion-
ization models but consistent with other AGN observa-
tions (Seaton 1978, Dumount et al. 1998).
34- C'CF analysis
An estimate of size and structure of the broad-line region
can be obtained from the cross-correlation function (CCF)
of a continuum light curve with emission line light curves.
We cross-correlated the 5100 A continuum light curve with
all our emission line light curves (Fig. 4) using an in-
terpolation cross-correlation function method (ICCF) de-
scribed by Gaskell & Peterson (1987). In Fig. 9 we plot
the cross-correlation functions of the individual emission
line light curves of HeIIA4686, HeIA5876, H/3 and Ha with
the continuum light curve. The cross-correlation functions
of HcIA4471 and HcIA5876 are identical within the errors;
therefore, only one curve is shown in the plot.
First of all we determined an error of the centroids of
the ICCFs by averaging the centroids T cent that were cal-
culated for fractions of the peak ranging from 35% to 90%
of the maximum value of the cross-correlation functions.
Then we estimated the influence of two principal sources
of cross-correlation uncertainties namely flux uncertainties
Table 5. Cross-Correlation Lags
Line
'Tcent
[days]
(1)
(2)
HeIIA4686
9.41,1;!
H/3
39.9t 3 |g
HeIA5876
59.6l3g;g
HeIA4471
62.7±f 7 5 7
Ha
81-6±I?;t
in individual measurements and uncertainties connected
to the sampling of the light curves. We used a method
similar to that described by Peterson et al. (1998b). We
added random noise to our measured flux values and cal-
culated the cross-correlation lags a large number of times.
Due to the large variability amplitudes of Mkn 110 these
uncertainties were of lower weight compared to those in-
troduced by the sampling of the light curves. The sampling
uncertainties were estimated by considering different sub-
sets of our light curves and repeating the cross-correlation
calculations. Typically we excluded 37% of our spectra
from the data set (cf. Peterson et al. 1998b). In Table 5
we list our final cross-correlation results together with the
total error.
Considering the entire observing period we got a lag
of the H/3 light curve of 39.9^g 3 5 2 days. Peterson et al.
(1998a) obtained for a similar extended observing cam-
paign a lag of 31.6l7;3 days. However, they claim that
their best lag estimate derived from an observing period
of 123 days yielding the smallest error is about 19.5tg'J
days.
K. Bischoff & W. Kollatschny: Strong optical line variability in Mkn 110
0.8
0.6
0.4
0.2
HoIIAfiftfi (1
neii4 , DOD i
Hol^R7fi f
neioo fo i
- H/S |
1
— Ha [
1?
1
jj
H
: J
W -
s />
f /si
\\ ' s -
f ///
V
Table 6. Mean and rms line widths (FWHM)
-6000
-3000
v [km/s]
3000
6000
Fig. 10. Mean profiles of HcIIA4686, HeIA5876, H/3 and
Ha.
0.8
0.6
0.4
0.2
1 1 1 1
i i i
- HeII4686
- HeI5876
H/3
1 h.
- — Ha
i in "i- (■■
J r\ V,
A \
-
i 1
: / !
S J 1
1
\V n j
\ :
- r JJ^V'/
j
-6000 -3000 3000 6000
v [km/s]
Fig. 11. The rms profiles of HeIIA4686, HeIA5876, H/3 and
Ha.
3.5. Line profiles and their variations
Normalized mean and rms profiles of HeIIA4686,
HcIA5876, H/3 and Ha lines are shown in Figs. 10 and
11. The rms profile is a measure of the variable part in
the line profile. There is a very broad line component in
the mean and rms profiles especially to be seen in the Hell
line. Even apart from this very broad component the mean
and rms profiles of the individual lines are different with
respect to their shape and full width at half maximum
(FWHM). In Table 6 we list the widths of the mean and
rms profiles. The mean and rms H/3 profiles are more sim-
ilar to the HeIA4471 profiles than to Ha. The rms profile
of Ha e.g. is significantly narrower than the profile of H/3.
The profiles of the HeIA4471 line are more noisy than the
other ones. They are identical to those of the HeIA5876
line within the errors.
Line
(1)
FWHM(mean)
[km s" 1 ]
(2)
FWHM(rms)
[km s^ 1 ]
(3)
HeIIA4686
2720 ± 100
4930 ± 200
H/3
1670 ± 50
2010 ± 100
HeIA5876
1640 ± 50
2200 ± 100
Ha
1580 ± 50
1540 ± 100
0.15
0.1
0.05 -
Ha JD7574
Ha JD9237
-6000
-3000
3000
6000
v [km/s]
Fig. 12. Ha spectra taken at the first minimum state of
February 1989 (solid line) as well as at the second mini-
mum state of August 1994 (dotted line) and their differ-
ence spectrum (dashed line).
All mean and rms profiles show a red asymmetry. The
asymmetry is mainly caused by a second line component
at v=1200 km s _1 . This second component does not vary
with the same amplitude as the main component. Further-
more, this second component was stronger during the first
half of our campaign from 1987 until January 1992 than
during the second half of the campaign. The Ha spec-
tra taken at the intensity minima of February 1989 and
August 1994 are plotted in Fig. 12. The additional com-
ponent centered at v=1200 km s _1 is clearly to be seen.
The mean spectra of the first half of our campaign are
broader by 400 - 500 km s" 1 (FWHM) than those of the
second half because of this component.
There is an independent very broad component present
in the mean and rms Hell profiles (Figs. 10 and 11). This
very broad component exists in addition to the broad com-
ponent. There is no transition component visible in the
profile. The peak of this very broad profile component is
redshifted by 400±100 km s _1 with respect to the narrow
lines. This shift was measured in the difference spectra
(cf. Fig. 13). This very broad component is the strongest
contributor to the Hell variability as can be seen from
the rms profile. The very broad component is visible in
the Balmer line profiles also, especially at high contin-
8
K. Bischoff & W. Kollatschny: Strong optical line variability in Mkn 110
-6000 -3000 3000
v [km/s]
Fig. 13. The very broad components of the Hell and H/3
profiles taken in January 1992 after subtraction of our
minimum profile taken in October 1988.
uum stages (see Fig. 3). The Hell and H/3 profiles taken
in January 1992 are shown in more detail in Fig. 13. We
subtracted the minimum profile taken in October 1988 to
remove the narrow line component. The Hell line intensity
has been divided by a factor of 1.3 for direct comparison
with the very broad H/3 profile. Immediately one can see
the striking similarity. The blue wing of the H/3 profile is
stronger than that of the Hell profile because of the blend-
ing with the red wing of the Hell line. The very broad line
component has a full width at zero intensity (FWZI) of
12 000 km s" 1 .
4. Discussion
Mkn 110 is one of the very few Seyfert galaxies with spec-
tral variability coverage over a time interval of ten years.
Different continuum ranges show different variability am-
plitudes; this holds for different optical emission lines, too.
But the mean fluxes of the continuum and of all emission
lines remain nearly constant integrated over time scales of
a few years (see Fig. 4) . There are considerable variations
over time scales of days to years. The strongest variability
amplitudes in the continuum shows the blue spectral range
(see Figs. 2 to 4). There are intensity variations of a fac-
tor of > 5. The strongest amplitudes in the blue spectral
range might be explained by a greater share of the non-
thermal continuum with respect to the underlying galaxy
continuum.
The optical line variations of H/3 are very strong in
comparison to other Seyfert galaxies (e.g. Peterson et al.
1998a). The HcIIA4686 line shows the strongest variations
of nearly a factor of 8 within two years. On the other
hand the H/3 and the continuum(A5100) vary only by a
factor of 1.7 and 3.0 respectively within the same time in-
terval. Apart from the variation of the HeIIA4686 line in
NGC 5548 in 1984 (Peterson & Ferland 1986) these are
the strongest optical line variations within such a time
interval. In the case of Mkn 110 we can show that the ap-
pearance of the very broad HeIIA4686 and H/3 component
(see Figs. 3 and 13) is not a unique event in the accretion
rate. It is connected to a very strong ionizing continuum
flux as can be seen from the light curves. The intensity ra-
tio HcIIA4686/H/3 comes to a value of 1.3 (see Fig. 13) in
the very broad line region. Such a line ratio is still in cor-
respondence with photoionization of broad emission-line
clouds in quasars (Korista et al. 1997).
The very broad line region (VBLR) originates close to
the central ionizing source at a distance of about 9 light
days. It is not connected to the "normal" BLR. As can
be seen from the line profiles there exists no continuous
transition region between these BLRs. The center of the
VBLR line profiles is shifted by 400±100 km s" 1 with
respect to the "normal" BLR profiles (Figs. 10, 11, 13).
Apart from this VBLR component we could show that
the line profiles of the Balmer and Hel lines arc similar
but not identical. The Ha line profile is narrower than
the H/3 profile. Besides the cross-correlation results this
is an independent indication that these two lines do not
originate in exactly the same region.
The observed Lya/H/3 ratio comes to a value of 11.0
in Mkn 110. This is about a factor of two higher than
the mean observed Lya/H/3 ratio in Seyfert galaxies (Wu
et al. 1983). Photoionization models of Kwan & Krolik
(1981) result in Lya/H/3 ~ 10 without the presence of
dust. Therefore, dust may not play an important role in
the BLR of Mkn 110. The Balmer decrement in Mkn 110
varies as a function of the ionizing continuum flux. This
might be explained by radiative transfer effects rather
than by variation of dust extinction.
The profiles of the broad emission lines in Mkn 110
are neither symmetric nor smooth (Figs. 10, 11). This is a
further indication that the broad-line regions in AGN are
structured as e.g. in NGC 4593 (Kollatschny & Dietrich
1997). In Fig. 12 it is shown that during the first half of
our campaign a red line component was present in the Ha
spectra. This component was not visible during the second
half of the campaign.
The size of the H/3 line emitting region (r = 40 Id
corresponding to 1.0 10 17 cm) and the optical continuum
luminosity is compared to those of other Seyfert galaxies.
The continuum luminosity amounts to L5100 = 4.4 10 39
erg s _1 A -1 . In this case we used H = 100 km s _1 Mpc -1
in order to compare directly the radius and luminosity of
Mkn 110 with those of other Seyfert galaxies compiled by
Carone et al. (1996). The values of Mkn 110 fit nicely into
the general radius-luminosity relationship for the broad-
line regions in Seyfert galaxies. The Balmer line emitting
region as well as the luminosity of Mkn 110 are arranged
in the upper region of the radius-luminosity plane close to
the galaxies Mkn 590 and Mkn 335.
K. Bischoff & W. Kollatschny: Strong optical line variability in Mkn 110
9
Table 7. Virial mass estimations
Line FWHM(rms) R M
[km s- 1 ] [Id] [10 7 M Q ]
(1) (2) (3) (4)
HeIIA4686
4930
±
200
9.4±
5.6
12.8
6-7l 4 9 ;?
HeIA5876
2200
±
100
59 6 +43 1
H/3
2010
±
100
39.91
33.2
9.5
4.7±*i
Ha
1540
±
100
8i.6l 2 3?;?
There is a trend that the broader emission lines origi-
nate closer to the center (see Table 7). A similar trend was
found for NGC 5548, too (Kollatschny & Dietrich 1996).
The central mass in Mkn 110 can be estimated from
the width of the broad emission line profiles (FWHM)
under the assumption that the gas dynamics are dom-
inated by the central massive object. Furthermore, one
needs the distance of the dominant emission line clouds
to the ionizing central source (e.g. Koratkar & Gaskcll
1991, Kollatschny & Dietrich 1997). We presume that the
characteristic velocity of the emission line region is given
by the FWHM of the rms profile and the characteristic
distance R is given by the centroid of the corresponding
cross-correlation function:
M = 'to^G^R.
In Table 7 we list our virial mass estimations of the central
massive object in Mkn 110. Altogether we determine a
central mass of:
M = 6.4lJ^10 7 M Q .
We can independently estimate an upper limit of the
central mass if we interpret the observed redshift of the
very broad Hell component (Az — 0.0013) as gravita-
tional redshift (e.g. Zheng & Sulentic 1990):
M = c 2 G~ 1 RAz
Again we presume that this line component originates at
a distance of 9 Id from the central ionizing source. We
derive an upper limit of the central mass of
M = 2.1 1O 8 M
This second independent method confirms the former
mass estimation.
5. Summary
Mkn 110 shows strong variations in the continuum and
in the line intensities on time scales of days to years. The
continuum - especially the blue range - varies by a factor of
3 to 5 on time scales of years. The Balmer line intensities
vary by a factor of 2.5 while the HcIIA4686 line shows
exceptionally strong variations by a factor of 8.
We cross-correlated the light curves of the emission
lines with those of the continuum. The emission lines orig-
inate at distances of 9 to 80 light days from the central
source as a function of ionization degree.
Not only the line intensities but also the line profiles
varied. We detected a very broad line region VBLR com-
ponent in the high intensity stages of the Balmer and
Hell lines. This region exists separated from the "normal"
broad-line region at a distance of only 9 light days from
the central ionizing source.
We derived the central mass in Mkn 110 using two
independent methods.
Mkn 110 is a prime target for further detailed variabil-
ity studies with respect to the line and continuum variabil-
ity amplitudes as well as with respect to the short-term
variations.
Acknowledgements. We thank M. Dietrich, D. Grupe, and U.
Thiele for taking spectra for us. We are grateful to M. Dietrich,
E. van Groningen, and I. Wanders who made available some
software to us. This work has been supported by DARA grant
50 OR94089 and DFG grant Ko 857/13.
References
Antonucci R.R.J. , Cohen R.D., 1983, ApJ 271, 564
Carone T.E., Peterson B.M., Bechtold J., et al., 1996, ApJ 471,
737
Dietrich M., Kollatschny W., Alloin D., et al., 1994, A&A 284,
33
Dumount A-M., Collin-Souffrin S., Nazarova L., 1998, A&A
331, 11
Gaskell CM., Peterson B.M., 1987, ApJS 65, 1
Hutchings J.B., Craven S.E., 1988, AJ 95, 677
Kollatschny W., Dietrich M., 1996, A&A 314, 43
Kollatschny W., Dietrich M., 1997, A&A 323, 5
Koratkar A., Gaskell M., 1991, ApJ 370, L61
Korista K, Baldwin J., Ferland G, Verner D., 1997, ApJS 108,
401
Kwan J., Krolik J.H., 1981, ApJ 250,478
Nelson C.H., MacKenty J.W., Simkin S.M., Griffith R.E., 1996,
ApJ 466, 713
Osterbrock D.E., 1989, Astrophysics of Gaseous Nebulae and
Active Galactic Nuclei, University Science Books, Mill Val-
ley
Peterson B.M., Ferland G.J., 1986, Nature 324, 345
Peterson B.M., Balonek T.J., Barker E.S., et al, 1991, ApJ
368, 119
Peterson B.M., Wanders I., Bertram R., et al., 1998a, ApJ 501,
82
Peterson B.M., Wanders I., Home K, et al., 1998b, PASP 110,
660
Petrosian A.R., Saakian K.A., Kachikian E.E., 1978, As-
trofizika 14, 69
Robinson A., 1994, The LAG Spectroscopic Monitoring Cam-
paign: An Overview. In: Gondhalekar P.M., Home K, Pe-
terson B.M. (eds.) Reverberation Mapping of the Broad-
Line Region in Active Galactic Nuclei, ASP Conf. Ser. 69,
San Francisco, p. 147
Rodri'guez-Pascual P.M., et al., 1997, ApJS 110, 9
10
K. Bischoff & W. Kollatschny: Strong optical line variability in Mkn 110
van Groningen E., Wanders I., 1992, PASP 104, 700
Weedman D.W., 1973, ApJ 183, 29
Wu C.-C, Boggess A., Gull T.R., 1983, ApJ 266, 28
Zheng W., Perez E., Grandi S.A., Penston M.V., 1995, AJ 109,
2355
Table 2. Continuum and integrated line fluxes
Julian Date 5100 A 4265 A 3750 A Ha H/3 HeIIA4686 HeIA4471 HeIA5876
2 440 000+
(1) (2) (3) (4) (5) (6) (7) (8) (9)
6828
3.14
±
0.09
3.97
±
0.1
9.0 ± 0.2
1656
±
24
477.9 ±
9.8
108.3 ± 5.2
30.8
±
5.2
68.7
±
5.4
7229
2.57
±
0.08
2.83
±
0.1
5.2 ± 0.2
1499
±
22
389.3 ±
8.9
58.1 ± 4.2
20.6
±
4.4
54.0
±
4.7
7438
1.77
±
0.07
1.96
±
0.1
2.9 ± 0.4
1112
±
18
287.3 ±
7.9
37.2 ± 3.7
17.6
±
4.2
33.1
±
3.7
7574
3.06
±
0.09
3.53
±
0.1
6.6 ± 0.4
1439
±
21
347.0 ±
8.5
60.5 ± 4.2
26.7
±
4.9
51.1
±
4.6
7663
4.43
±
0.12
5.27
±
0.1
12.5 ± 0.4
1744
±
24
475.1 ±
9.8
159.2 ± 6.2
25.8
±
4.8
61.8
±
5.1
7828
1.98
±
0.07
1.99
±
0.1
1520
±
22
357.9 ±
8.6
30.6 ± 3.6
22.5
±
4.6
41.7
±
4.1
8632
5.93
±
0.15
7.56
±
0.1
15.1 ± 0.4
1995
±
27
624.5 ±
11.2
237.6 ± 7.8
41.7
±
5.9
103.
±
7.2
8812
3.24
±
0.09
2091
±
28
494.9
±
9.9
101.2 ± 5.0
36.1
±
5.5
95.6
±
6.8
8860
2.86
±
0.09
1673
±
24
421.4
±
9.2
80.7 ± 4.6
35.4
±
5.5
70.7
±
5.5
8862
3.26
±
0.10
1782
±
25
454.9
±
9.5
92.1 ± 4.8
32.6
±
5.3
81.7
±
6.1
8864
3.30
±
0.10
3.00
±
0.2
4.8 ± 0.2
1703
±
24
432.1
±
9.3
80.7 ± 4.6
32.0
±
5.2
79.7
±
6.0
9078
3.15
±
0.09
1692
±
24
468.5
±
9.7
91.2 ± 4.8
30.9
±
5.2
80.4
±
6.0
9080
3.12
±
0.09
1702
±
24
453.9
±
9.5
88.2 ± 4.8
23.9
±
4.7
80.4
±
6.0
9083
3.22
±
0.09
1764
±
25
467.7
±
9.7
86.6 ± 4.7
27.8
±
4.9
77.6
±
5.9
9123
3.47
±
0.10
1678
±
24
482.3
±
9.8
1 59 4 + 6 2
40.9
±
5.9
87.3
±
6.4
9237
1.82
±
0.07
2.01
±
0.2
4.1 ± 0.4
1197
±
19
310.2
±
8.1
65.5 ± 4.3
25.6
±
4.8
59.0
±
4.9
9419
1.91
±
0.07
288.0
±
7.9
81.9 ± 4.6
9595
1.53
±
0.06
876 ± 16
266.8
±
7.7
46.5 ± 3.9
12.0
±
3.8
38.4
±
3.9
9776
4.00
±
0.11
4.95
±
0.1
1432 ± 21
414.6
±
9.1
133.9 ± 5.7
29.1
±
5.0
68.1
±
5.4
9785
4.17
±
0.11
4.57
±
0.1
1430 ± 21
411.8
±
9.1
135.9 ± 5.7
27.9
±
5.0
72.5
±
5.6
9786
3.95
±
0.11
4.30
±
0.1
1600 ± 23
412.9
±
9.1
128.4 ± 5.6
22.6
±
4.6
74.9
±
5.7
9787
3.96
±
0.11
4.44
±
0.1
1529 ± 22
433.6
±
9.3
128.2 ± 5.6
30.9
±
5.2
72.1
±
5.6
9811
3.62
±
0.10
4.36
±
0.1
1485 ± 22
441.1
±
9.4
95.1 ± 4.9
35.0
±
5.4
69.3
±
5.5
9870
4.86
±
0.13
1558 ± 23
445.9
±
9.5
174.2 ± 6.5
75.8
±
5.8
Continuum fluxes (2) - (4) in 10 15 erg sec 1 cm 2 A
Line fluxes (5) - (9) in 10 -15 erg sec -1 cm~ 2 .